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Abstract. The static transverse and longitudinal 
correlation functions (CF) of a 3-dimensional ferromag- 
net are calculated for the exactly solvable anisotropic 
spherical model (ASM) determined as the limit D — + oo 
of the classical D-component vector model. The re- 
sults are nonequivalent to those for the standard spher- 
ical model of Berlin and Kac even in the isotropic case. 
Whereas the transverse CF has the usual Ornstein- 
Zernike form for small wave vectors, the longitudinal CF 
shows a nontrivial behavior in the ordered region caused 
by spin-wave fluctuations. In particular, in the isotropic 
case below T c one has S zz (\sl) oc \/k (the result of the 
spin- wave theory) for k < K m oc T c — T. 

PACS: 75.10.Hk; 75.40.Cx 



The spherical model (SM) of Berlin and Kac was be- 
lieved many years to be the only model in the statistical 
theory of magnetism, which is exactly solvable in 3 di- 
mensions. The other model, which is the limit D — > oo 
of the classical D-component vector model introduced by 
Stanley |Q| , has the same partition function in the homo- 
geneous case 0], and for this reason the latter was not 
considered as an independent model by many researchers 
and received relatively little attention. 

However, the D = oo model is nonequivalent to the 
standard spherical model of Berlin and Kac and goes be- 
yond it in many respects. First, it can be easily gener- 
alized for anisotropic systems Q. Second, it does not 
use the global spin constraint, which leads to unphys- 
ical results in spatially inhomogeneous situations. In 



particular, in the SM the Curie temperature T c of a 4- 
dimensional ferromagnetic film with free boundary con- 
ditions, which is infinite in 3 dimensions and finite in the 
4th dimension, has been found by Barber and Fisher || 
to be a non-monotonous function of the number of its lay- 
ers N . To the contrast, an improved version of the spher- 
ical model using separate spin constraints in each layer || 
leads in the case mentioned above to the monotonically 
increasing T C (N), as it should be. However, even this 
improved version of the SM fails on physically meaning- 
ful 3-dimensional ferromagnetic films, since the latter are 
2-dimensional objects and without anisotropy their T c is 
zero for any finite thickness. An adequate description of 
spatially inhomogeneous and low-dimensional ferromag- 
netic systems on the "spherical" level can be archieved 
only with the use of the anisotropic D — oo model, which 
can be called also anisotropic spherical model (ASM). 
This model was recently applied to study the dimensional 
crossover of T C (N) of anisotropic ferromagnetic films Q| 
and the effects of thermal fluctuations in the Bloch-wall 
phase transition || . 

Here it will be shown that the ASM deviates from the 
standard SM even in the spatially homogeneous case and 
even in the isotropic limit, if the spin-spin correlation 
functions are concerned. Indeed, below T c in the D = oo 
model there are two different — longitudinal and trans- 
verse — correlation functions, whereas there is only one 
CF in the SM. A less trivial reason for the difference be- 
tween two models relies on the fact that a wave-vector- 
dependent CF is proportional to the appropriate suscep- 
tibility, which is the linear response to a spatially inho- 
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mogeneous sinusoidal field. In this case, as was argued 
above, the global spin constraint used by the SM modifies 
the results. As we shall see, the k-dependent longitudi- 
nal CF of the ASM has a nontrivial non-Ornstein-Zernike 
form below T c , which is similar to that following from the 
lowest-order spin- wave theory for T -C T c . For 

the uniaxially anisotropic ferromagnetic model the ASM 
yields finite longitudinal and transverse susceptibilities 
and correlation lengths in the ordered region. 

The anisotropic generalization of the classical D- 
component vector model of Stanley pj can be described 
by the Hamiltonian 
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where |m, | = 1 and r\ a < 1 are anisotropy coefficients. 
Here we consider the uniaxial model with 771 = 77^ = 1 
and rj a = r\ < 1 for a > 2. The model (]l|) can be con- 
veniently treated by the classical spin diagram technique 
p2| , ^3|Jl4| , which allows classification of diagrams in pow- 
ers of l/D for D > 1 Hi . The equation of state m(H, T) 
of the anisotropic spherical model, where to = (m z ), 
is contained in the diagram series corresponding to the 
self-consistent Gaussian approximation (SCGA) [fl2Lfl4j, 
which becomes exact in the limit D — ► 00. In terms of 
the dimensionless variables 9 = T/T^ FA and h = H/Jo, 
where T C MFA = Jo/D is the mean- field transition tem- 
perature and Jo is the zero Fourier component of the 
exchange interaction, one comes to the system of equa- 
tions for the magnetization m and the gap parameter G 

BR 



G 



m + h ' 
6GP(r)G) = 1 - to 2 . 



Here 



P(X) = v 



(2tt) 3 1 - AA q 
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(3) 



where vq is the unit cell volume and A q = J q / Jo- In the 
long- wavelength limit A q = 1 — aq 2 , where a ~ a 2 , and 
Oo is the lattice spacing. For 3-dimcnsional lattices with 
the nearest neighbour interactions the integral P(X) has 
the following properties: 



P(X) 



l + X 2 /z, 
W-c {l-Xf/ 2 , 



A< 1 
1-X«1, 



(4) 



where z is the number of nearest neighbors, W is the 
Watson integral, and Co = 7Jo/(47ra 3 / 2 ). For the simple 



cubic (sc) lattice vo 



— ^3 



and a = a 2 ,/ 6, hence Co 



(2/7r)(3/2) 3 / 2 . The solution of the system of equations 
(^|) simplifies for zero field, where below T c one has G = 1 
(the zero spin-wave gap) and 
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It can be seen that in the isotropic case, ry = 1, the 
value of the phase transition temperature 6 C reduces to 
the well-known result 6 C = 1/W @. Using (§) one can 
calculate the longitudinal susceptibility Xz = dm/dH. 
The zero-field reduced susceptibility, \ z = JqXz, has the 
form 
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where P'(X) = dP(X) / dX and G satisfies the equation 
0GP(rjG) — 1 above 8 C . Solving this equation near 9 C in 
the linear approximation in 1 — G <C 1, using (|5|), and 
introducing the reduced temperature variable 
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one can rewrite (0) in the form 
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i] < 1 
1 - f] < 1. 
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The first line of (|J) is valid for the weakly-anisotropic 
model, 1 — rj <C 1, in the narrow temperature interval 
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whereas the second one is valid in the whole region be- 
low 9 C . The latter is the reason to define e in the non- 
standard form (Q). It can be seen from (H) that the crit- 
ical behaviour of the longitudinal susceptibility in the 
ASM is the same as that in the mean field approxima- 
tion (MFA), including the famous ratio of critical am- 
plitudes 2. This equivalence holds for all the critical in- 
dices |^2[ , since for 77 < 1 the square-root singularity of 
P{X) at X = 1 is suppressed. In the extreme case of the 
"spherical Ising model", i] = 0, the mean field approx- 
imation for the ASM becomes exact, since the fluctua- 
tions of the transverse spin components die out with the 
transverse coupling 77 and the influence of longitudinal 
fluctuations vanishes in the spherical limit, D — > 00. For 
the anisotropic model, rj < 1, the quantity \ z is finite in 
the ordered region, vanishing at 6 = and diverging at 
9 = 6 C . On the contrary, the longitudinal susceptibility 
of the isotropic model, 77 = 1, diverges as \z (— e)~ 2 
above 9 C and is infinite in the whole region below 9 C due to 
the Goldstone-mode fluctuations. In this case for h 1 
one has 
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FIG. 1. Temperature dependences of the longitudinal and 
transverse susceptibilities of the uniaxial spherical model. 



(11) 



which coincides with the result of the standard spher- 
ical model and is in accord with the prediction of the 
lowest-order spin- wave theory |l5| for 9 <C C . A simi- 
lar expression for a general D-component vector model, 
which contains the additional factor 1 — 1/ D, was derived 
by Fisher and Privman [T^ | from the scaling arguments. 
Measuring the square-root singularity of magnetization 
( pd| ) in systems with spontaneous breaking of a contin- 
uous symmetry is a difficult task, and it was done only 
resently on EuS [|l7j and on EuO Q . 

In the following we shall study the correlation functions 



S aa (r) = (Am Q (0)Am Q (r)>, 



« = 1,2, 



,A (12) 



which are translationally invariant in the spatially homo- 
geneous case considered here. The fluctuation-dissipa- 
tion theorem relates wave-vector-dependent CFs to the 
appropriate susceptibilities as XaO*-) — /3S' Q , Q (k), where 
f3 = 1/T. Both longitudinal and transverse CFs enter 
the self-consistent Gaussian appro ximation and were (ap- 
proximately) calculated in [jl2| ^ p~4|| . Since the SCGA be- 
comes exact in the limit D — > oo, the transverse CFs 
describing in the SCGA fluctuations of D — 1 transverse 
spin components become exact, too, and with the use of 
the results of ifhl one obtains 



*±(k) = £sj.(k) = G 



(13) 



in the whole temperature range [cf. (§]. One can see 
that for h = below 8 C , where G = 1, the homoge- 
neous transverse susceptibility, x± = x_l(0), diverges in 
the isotropic limit, r] — > 1, and remains a finite constant 
for the uniaxially-anisotropic model. Comparing (ra) and 
(|l3|), one finds that both susceptibilities become equal 
to each other at the characteristic temperature slightly 
below C , which is determined by e = e* [see (0) and 
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FIG. 2. (a) The compact part A zz (k) of the longitudinal 
CF 5 zz (k); (b) block summation for one-site spin cumulants 
A... renormalized by Gaussian fluctuations of the molecular 
field; (c) ladder equation for the transverse four-spin correla- 
tion line Vk; (d) Dyson equation for the renormalized trans- 
verse interaction 7?/3J q . Unlabeled small circles denote trans- 
verse spin components. 



(|Io|)]. The temperature dependences of the longitudinal 
and transverse susceptibilities for the weakly-anisotropic 
spherical model are represented in Fig. [|. One can see 
that even a small anisotropy, as 1 — T) = 1CU 2 , has a 
profound influence on Xz below 9 C . 

Unlike the transverse CF given by (|l3|), the exact lon- 
gitudinal CF, S zz (k), cannot be determined from the 
SCGA in the spherical limit, since for D — > oo the fluctu- 
ations of the single longitudinal spin component become 
nonessential and S zz (Js) disappears from the SCGA equa- 
tions. Therefore, the wave-vector-dependent longitudinal 
CF in the anisotropic spherical model should be consid- 
ered separately, which is the main purpose of this paper. 
With the help of the classical spin diagram technique 
©HUH S zz (k) determined by @ can be represented 
as 



S zz (k) 



A zz (k) 



l-A zz (k)/3J k 



(14) 



where A zz (k) is the compact (irreducible) part of 5 zz (k) 
given by the diagrams, which cannot be cut by the one 
longitudinal interaction line /3Jk- The quantity A zz (k) 
is in turn given in the limit D — ► oo by the set of dia- 
grams represented in Fig. p[ Such a choice of diagrams 
is based on the arguments of |l3| . More technical details 
can be found in [|19[ , where the transverse CF, S' QQ (k), 
a > 2, was calculated up to the first order in 1/D for 
low-dimensional ferro- and antiferromagnets in magnetic 
field. The analytical form of A zz (k) of Fig. ^a reads 



A zz (k) = A zz + A 2 QZ Vk. 
Here in the limit D — > oo one has 
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A-zz — A QQ + A QQZ — Aaajg^, (16) 

where a ^ (3 7^ z, 



G satisfies the system of equations (0) , and 



(17) 



D Dm 
£ = + m J ) = -(ft + m) = — (18) 

is the temperature-normalized molecular field. The quan- 
tity Vk in ( Jl5| ) is the solution of the ladder equation Fig. 
||c and has the form 
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where A QQ(3/3 is given by ([H]), 
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and the renormalized transverse interaction ry/JJq deter- 
mined by Fig. 0d reads 



?7/3 J q = 



ri/3 J q 



1 - A aa r)(3J q 



(21) 



The factor Z) — 1 in ( pp[ ) results from the summation over 
transverse spin components in Fig. p|a or in Eq. ( |l5| ) , as 
well as in the ladder equation Fig. |^c! Such a factor does 
not appear, if one tries to take into account the similar 
diagrams for the transverse CF <Sj_(k) = 5 QQ .(k), a > 2, 
and thus these diagrams vanish in the spherical limit and 
-Sj_(k) has the trivial Ornstcin-Zcrnike form (p"3|). 

Combining now Eqs. ([l4|)-(pl|), one comes after sim- 
plifications to the final result 
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(23) 



The result similar to (22), in which, however, only the 
last term is present, was obtained earlier ||Jl(],[llJ within 
the lowest-order spin- wave theory well below T c . On the 
other hand, for the standard spherical model the cor- 
relation function has the trivial Ornstein-Zernike for all 
temperatures Pl^oft. The integral in (|23| ) can be easily 
calculated for k = and in the corner of the Brillouin 
zone, k = b, for the sc lattice, where Ab-q = — A q . The 
result has the form 



r k = P(r)G) 



I(r)G), k = 



(24) 



where I(X) is given by (0). One can see that xj 1 (0) of 
( 22) is in accord with the previously obtained expression 
(W for h = below 9 C , where G = 1. Under the same 
conditions xJ 1 (b) = 8/(29 c ), which can be used to con- 
trol numerical calculations. For the weakly anisotropic 
model in the case 1 — r\G <C 1 the integral ( |23] ) can be 
calculated analytically for small wave vectors, a^k <C 1, 
which results in 
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(25) 



The last term in (g2j) modifies the k-dependence of Xz (k) 
below 9 C in the gapless case nG = 1 in the range of small 
wave vectors, which shrinks if 9 C is approached from be- 
low. Instead of the Ornstein-Zernike form X*(k) oc l/k 2 
one has (cf. |M0) 



Xz(k) OC -, k < K r , 
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(26) 



The latter defines a length scale £ m , which exists only in 
the ordered region (m > 0) and diverges at 9 C due to the 
vanishing of magnetization (0): 
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(27) 



and hence £„ 



38 „ 



For the sc lattice coce 1 / 2 
The length £ m is analogous to the "bare", i.e., the mean- 
field correlation length below T c , which follows, in partic- 
ular, from the Landau- Ginzburg phenomcnological free 
energy. This analogy is, however, not complete, since £ m 
diverges at the actual transition temperature T c and not 
at T C MFA (6> = 1). The crossover of x*(k) at k ~ K m for 
temperatures slightly below T c was described earlier with 
the help of the renormalizati on g roup approach pi]] . The 
k-dependences of xj 1 (k) of ( ^2|) in the isotropic case ob- 
tained with the help of numerical integration in (^3|) are 
represented in Fig. |§| for different temperatures below 
T c . Curves of such a type for Heisenberg model systems 
such as EuO and EuS could be observed, in principle, 
in neutron scattering experiments, but such experiments 
were not carried out up to now. 

Our next task is to calculate explicitly the real-space 
correlation functions ( ]l2| ) in the small-anisotropy case 
1 — f] <C 1. Since the CFs themselves are proportional 
to 1/D and thus vanish in the spherical limit [see, e.g., 
dl3|)], it is more convenient to deal with the appropriate 
susceptibilities, which are given by 



(28) 



In the transverse case using (|l^) one comes to the well- 
known result 



t V 1 — f]G 



(29) 
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FIG. 3. The wave vector dependences of the inverse longi- 
tudinal susceptibility Xz" 1 (k) of the isotropic spherical ferro- 
magnet along the [111] direction of k for different tempera- 
tures in the ordered region. 

for r > oo. One can see that in the isotropic case, r\ = 1, 
the transverse correlation length, £ c j_, if infinite in the 
whole region below 9 C , where G = 1. In more complicated 
situations, as for the longitudinal susceptibility below 9 C , 
(p2|), the correlation length can be determined as £ cz = 
1/k c , where n c is the singularity point of Xz(k) on the 
imaginary axis k = in, which is closest to the origin. This 
leads with the use of ( p5| ) to the transcedental equation 
for u c = £ c ±/(2£ cz ) having the form 



u c In ■ 



1 



1 



Cm 



2m 2 



(30) 



0c o -y/l - V 

where £ m is given by (|27]). This equation coincides with 
Eq. (4.19) of ||, which determines the temperature- 
dependent width Sl — 2£ cz of the (linear) domain wall in 
the uniaxial spherical model. The asymptotic solutions 
of (|30|) read 



y/a/2, a < 1 (slightly below 9 C ) 
l-2e~ a , a>l (far below 9 C ). 



(31) 



Note that the transition between two regimes in (|3l| ) oc- 

The explicit results 



curs at e ~ e* ~ yl - r\ [see ( |l 0[ ) 
for the correlation length £ cz itself, including that above 
9 C , read 
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sm 



< -e< e* 

< e < e* ( 32 ) 



where I (if) is given by the second line of (j^. The numer- 
ically calculated temperature dependences of the longi- 
tudinal and transverse correlation lengths are illustrated 
in Fig. [| 
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FIG. 4. Temperature dependences of the transverse, £ c ±, 
and longitudinal, £ cz , correlation lengths, as well as the char- 
acteristic length £ m , in the uniaxial spherical model. 



To calculate the Fourier-transform \z (r) below 9 C with 
the help of (pgj), (p2[), and (p5|), it is convenient to per- 
form at first the integration over the angle variables and 
then to deform the ^-integration contour in the complex 
plane. In this way one comes to the well-behaved expres- 
sion 



2ac a 1 / 2 



du- 



-r/tc 



l-ui 



a(l-u 2 c ) + 2u 2 c 



wln^i+i 



(33) 



where the first term is the contribution of the pole deter- 
mining the longitudinal correlation length ^ cz [see (|3C|)] 
and the second one is the integral along the cut of the 
arctan- function of (p5|). The expression above simplifies 
in various limiting cases. For a <C 1 (e <C e*), the cut 
term of (33) is small as a in comparison to the pole term 
and, additionally, it decays as exp(— 2r/£ c jJ, i.e., much 
faster than the first term of (|33|). Thus, in this case only 
the pole term in ( p3| ) is relevant, and with the use of ( |3l| ) 
one obtains 



1 11 

Xz(r) = e r/4 ", 



(34) 



where ^ cz is given by the the middle line of (|3^). In the 
opposite case, a ^> 1 (e ^> e*), the pole sticks to the 
beginning of the cut, u = 1, and the amplitude of the 
pole term becomes exponentially small because of the 
factor 1 — u 2 [see ©]. Thus, in this case only the cut 
term in ( |33| ) is essential, and the result can be simplified 
to 



Xz(t) = 



2c n a 1 / 2 



-r/tc 



dx 



■ exp 



2r 



(35) 
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where £ cz is given by the the lower line of (p2[). 
asymptotic forms of (p5|) are (B4h for r <g; £ m and 



2m 2 /£ 



Mr) = -^r - 



The 



(36) 



for r 3> £ m . One can say that for the isotropic model, 
where the true correlation length £ cz is infinite below 8 C , 
the "bare" one, £ m of (p7|), still plays its role in some 
sence: For r > £ m the slow decay \ z (r) oc 1/r changes to 
a faster one, x z (r) °c 1/r 2 . 

As a conclusion, the anisotropic spherical model 
(ASM) considered here is a good exactly solvable "toy" 
model for classical spin systems, which can be success- 
fully applied in many situations where the standard 
spherical model of Berlin and Kac fails. More realistic 
models of phase transitions posess, naturally, the other, 
nonspherical, values of the critical indices, but this is, 
however, only a quantitative effect, which is less impor- 
tant in comparison to the profound role played by the 
gapless spin waves in the isotropic model below T c . The 
latter are properly taken into account in the ASM, which 
is thereby a very important step beyond the mean field 
approximation. 
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